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We study the heat conduction of the spin-1/2 XXZ chain in finite magnetic fields where mag- 
netothermal effects arise. Due to the integrability of this model, all transport coefficients diverge, 
signaled by finite Drude weights. Using exact diagonalization and mean-field theory, we analyze the 
temperature and field dependence of the thermal Drude weight for various exchange anisotropics 
under the condition of zero magnetization-current flow. First, we find a strong magnetic field depen- 
dence of the Drude weight, including a suppression of its magnitude with increasing field strength 
and a non-monotonic field-dependence of the peak position. Second, for small exchange anisotropics 
and magnetic fields in the massless as well as in the fully polarized regime the mean-field approach 
is in excellent agreement with the exact diagonalization data. Third, at the field-induced quantum 
critical line between the para- and ferromagnetic region we propose a universal low-temperature 
behavior of the thermal Drude weight. 



I. INTRODUCTION 



Transport properties of one-dimensional spin-1/2 
systems are currently at the focus of active re- 
search. This has been motivated by the exper- 
imental manifestation of significant contributions to 
the thermal conductivity originating from magnetic 
excitation9iiSiii^, stimulating intensive theoretical 
worliW,i"i",i2ii3,i4,i5,i6,i7,i8,i9_ ^^^^^^ theoretical 

effortsi'S'^iifliiiiiiii havc been devoted to the question of 
possible ballistic thermal transport in generic spin mod- 
els such as spin ladders, frustrated chains, and dimerized 
chains. Such ballistic transport would be characterized 
by a finite thermal Drude weight. Recent numerical and 
analytical studies indicate that in pure but nonintegrable 
spin models, the thermal Drude weight scales to zero in 
the thermodynamic limit implying that the thermal cur- 
rent is likely to have a finite intrinsic life-timeSiSiiSiiii^. 
In addition, the effects of extrinsic magnon scattering by 
phonons and/or impurities have been addressed in several 
works^&^^i^. For the integrablc XXZ model, the energy 
current operator is a conserved quantitySSiSi, leading to 
a finite thermal Drude weight. Its temperature depen- 
dence has been studied with exact diagonalizationZiSiS 
and Bethe ansatz techniquesiSiii and is well understood 
for arbitrary values of the exchange anisotropy at zero 
magnetic field. In this paper, we address the issue of 
thermal transport in the XXZ model in the presence of 
a finite magnetic field h. In this case, magnetothermal ef- 
fects become important and must be accounted for. The 
magnetothermal response itself has been studied by Louis 
and Gros in the limit of small magnetic fieldsiSi and re- 
cently also by Sakai and Kliimper in the low-temperature 
limiliiS. Here, we consider magnetic fields of arbitrary 
strength and we discuss the temperature dependence of 
the thermal Drude weight under the condition of zero 
spin-current flow. 



The Hamiltonian of the XXZ model reads 

^ = -^E I 2(^'^^'+i + H.c.) + A5f 5f+i - hSt^ (1) 

where TV is the number of sites, S^'^ are spin-1/2 op- 
erators acting on site I, and A denotes the exchange 
anisotropy. The exchange coupling J is set to unity in 
our numerical calculations. We focus on A > and pe- 
riodic boundary conditions arc imposed. 

The quantum phases and the spectrum of ^ are well 
understood, both as a function of exchange anisotropy A 
and magnetic field h. The reader is referred to Ref . 122 for 
a detailed summary and further references. Here we only 
repeat the main points. At zero magnetic field, the spec- 
trum of the Hamiltonian Eq. is gapless for |A| < 1 
and gapped for |A| > 1. The situation at finite mag- 
netic fields is summarized in the first four columns of 
Table |l| Three different cases are found: (i) the ferro- 
magnetic gapped state for h > hc2 = 1 + A (FM); (ii) 
the gapless or massless phase for h < hc2 = 1 -|- A and 
h > hci, and (iii) the antiferromagnetic, gapped state for 
A > 1 and h < hd (AFM). The hue h = hd starts at 
the SU{2) symmetric point A = 1, ft, = and hd grows 
exponentially slowly in the region A > 1, /i > 0. 







h 


mo 


r/j< 1 


(i) 


FM, gap 


h > hc2 


1/2 


A'th oc T^/^ exp(-G/T) 




Saturation 


h = hc2 


1/2 


Kth = const r^/^ 






hc2 = l + A 






(ii) 


Massless 


hci < h < hc2 




Ath cxT 


(iii) 


AFM, gap 


h < hci 








TABLE I: Magnetic phases of the XXZ model (see, e.g., 
Ref. 12^) and leading term of the thermal Drude weight Ath 
at low temperatures, mo is the average local magnetization 
at T = 0. G = G{h) denotes the gap in either the polarized 
state (i) or the massive antiferromagnetic regime (iii). In the 
polarized state (i), G{h)/J = h — hc2- 
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The fifth column of Table H] is a first account of our 
main findings for the low-temperature behavior of the 
thermal Drude weight, denoted by Kth in this paper. 
These results are now briefly summarized. One can ex- 
pect qualitative changes in the low-temperature behav- 
ior of the thermal Drude weight as the transition lines 
h = hci and h = hc2 are crossed. In particular, we focus 
on the transition from the gapless phase to the ferromag- 
netic state. In Sec. IIIII we will argue that for T/ J ^ 1, 
first, isTth cx T^^^ exp{—G/T) in the ferromagnetic state, 
G being the gap, and T temperature; second, Kth T 
in the massless phase; and third, Kth « T'^/^ along the 
line h = hc2- 

Regarding the antiferromagnetic state, there is cer- 
tainly also an exponentially suppressed Drude weight; 
see for instance Refs. andllTl for h = 0. However, the 
low-temperature region in this case and for h = hd is 
difficult to reach with the methods of the present paper. 
For a discussion of the low-temperature limit at vanish- 
ing magnetic field, we refer the reader to Refs. I^I^ITil and 
llTl Apart from the low-temperature behavior, this paper 
studies the field dependence of the thermal Drude weight 
in the phases (i) and (ii) at finite temperatures. 

The plan of this paper is the following. First, we dis- 
cuss the expressions for the transport coefficients and the 
current operators in Sec^l Second, in Sec. IIIII we per- 
form an analysis of the transport coefficients based on 
a Jordan- Wigner mapping of the spin system onto spin- 
less fermions. In this case, interactions at A 7^ will 
be treated by a Hartree-Fock approximation. Third, we 
present our results from exact diagonalization for A > 
in Sec. IIVI and compare them to the results from the 
Jordan- Wigner approach. The field and temperature de- 
pendence of the thermal Drude weight is discussed with 
a particular focus on the case of the Heisenberg chain. A 
summary and conclusions are given in Sec. |3 

II. TRANSPORT COEFFICIENTS 

Within linear response theory, the thermal and the 
spin current are related to the gradients V/i and VT of 
the field h and the temperature T hyM, 

\J2 J \ L21 L22 J \ -VT J 

where Ji = (ji) is either the thermodynamic expectation 
value of the spin current ji or the thermal current opera- 
tor j2, respectively. Lij denote the transport coefficients. 
At finite frequencies w, the coefficients Lij{uj) depend on 
the time-dependent current-current correlation functions 
via^^ 

or pea /■/? 

i..H = ^ X die*(-+'° )*y^ dr{jaj[t + zT)). (3) 

In this equation, r = for j = 1 and r = 1 for 
j = 2. /3 = 1/T is the inverse temperature and (.) de- 
notes the thermodynamic expectation value. Note that 



L12 = L21/T due to Onsager's relationS^. The real part 
of Lij (uj) can be decomposed into a 5-function at w = 
with weight Dij and a regular part Lj^^(w): 

ReL,y(c^) = Aj<5M + L^fM. (4) 

This equation defines the Drude weights Dij, for which 
a spectral representation can be given^ 

Aj (K T)= ^ ^ pn{n\ji\m){'m\jj\n) . (5) 

Ejn = En 

Here, p„ = exp(— /3£'„)/Z is the Boltzmann weight and 
Z denotes the partition function. In the exponent, r has 
to be chosen in the same way as in Eq. (Q. 

Let us now introduce the appropriate definitions of the 
current operators. The local current operators ji^i and 
j2,i satisfy the continuity equations 

jj.i+i - jj-.i = -^[H, rfj,;]; J = 1. 2 (6) 

where di i = S'f is the local magnetization density and 
d2.i = hi is the local energy density, respectively, with 
H ~ Tlii^i- At zero magnetic field, the total currents 
.7th[s] = E/ Jth[s],i are given bjiSLS^ 

N 

js = ijY.{S+Sr^,-S+^,Sr) (7) 
1=1 

jth = J'^S',-(S'/+i X J'j+a) (8) 
1=1 

with the definition 5"; = {Sf , Sf , ASi)to achieve a com- 
pact representation, while Si is defined as usual. Note 
that subscripts in brackets [.] refer to spin transport. 

At finite magnetic field, the proper set of current op- 
erators isi^ 

ji = is; j2 = jth - hjs. (9) 

Now, the crucial point is that, while the spin current js 
is only conserved in the XX case (A = 0), the current jth 
is conserved for all fields h and values of A, i.e., [H, jth] = 
(Refs. I20I and I21I) . Thus, it immediately follows from 
Eqs. Q and (jSJ that the Drude weights D12, D21, and 
D22 are finite for arbitrary fields h. 

Furthermore, one can show that the spin Drude weight 
Dii is also finite in the thermodynamic limit for /i ^ 0. 
We briefiy outline the proof along the lines of Ref. |^ 
Given a set of all conserved observables {Qi}, the spin 
Drude weight Dn can be written as 

Dii{h,T) = ^inlPj,) (10) 

where V is the projection operator in the Liouville space 
on the subspace spanned by all conserved quantities 
{Qi}- The brackets (.|.) denote Mori's scalar product; 
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see, e.g., Rcf. |26| for details. Restricting to a subset 
{Qm} C {Qi}, one obtains an inequalitjiSiiSi 



Dn{h,T)> 



TN 



E 



(11) 



providing a lower bound for the Drude weight Dnih, T). 
In the literature, this relation is often referred to as 
Mazur's inequalit3i2iiSi. Several authors^LS^ have used 
Eq. Hll|l to infer a finite spin Drude weight for the Heisen- 
berg chain, assuming broken particle-hole symmetry, or 
the presence of a finite magnetic field, respectively. More 
explicitly, only one conserved quantity is often considered 
in Eq. namely Qi ~ jth, which has a finite overlap 
(jibth) > with the spin current for h ^ 0. This finally 
proves Dii{h,T) > for /i 7^ 0. 

The main focus of this paper is on the case of purely 
thermal transport with a vanishing spin current, i.e., 
Ji = 0. We therefore arrive at a thermal conductivity 
K which is described by 

Re k{uj) = Kth{h, T)6{uj) + K,,g{uj) (12) 

where Kth{h, T) in terms of the Drude weights Dij reads 



Kti,{h,T) ^ D22{h,T) ~ f3 



Dn{h,T) 



(13) 



Exactly the same result for Kth{h,T) is obtained if a 
different choice of current operators and corresponding 
forces is made, e.g., jg and jth from Eqs. Q) and © 
(see Ref. 123). The expression for Kth{h,T), being fully 
equivalent to Eq. l(T^ . is then given by 



i^th(^T) = Ah(^T)-/3 



Ds{h,T) 



(14) 



Note that for h = 0, Kthih = 0, T) = Dth{h = 0,T). 
Therefore, two competing terms contribute to Kth{h,T) 
in Eq. (|14|l : the "pure" thermal Drude weight Dth and 
the " magnetothermal correction", PD^y^ ^/ Ds- Note that 
the magnetothermal correction might be suppressed by 
external scattering or spin-orbit coupling, breaking the 
conservation of the total magnetization of the spin sys- 
tem fRef. IT^. This is an open issue which may depend 
crucially on the particular material investigated in exper- 
imental transport studies. 

Let us now give spectral representations for t he quanti- 
ties D,{h,T), Dt^{h,T), and Ah,s(/i,T) (Refs. 



Lie quanti- 
l2ll29l3(TD 



Di\h,T) 

Ah,s(/i,r) 



7r/3' 



2[11 



TV 



E Pn\{m\3t\i[s\\n)\'^ , (15) 



TV 



7r/3 

Iv 



,(16) 
(17) 



The operator f = (1/2) I]i('S'i^'5'i+i +H-C-) is the kinetic 
energy. In the Eqs. IjlSfl . (|16|l . and 117|l . the magnetic 
field only enters via the Boltzmann weights p,, . The two 
expressions Dl and Dl^ are equivalent in the thermody- 
namic limit, but exhibit differences at low temperatures 
for finite system sizes^L^iiSiSi^. In this context, note 
that — D\ is the so-called Meissner fraction, which 
measures the superfluid density in the thermodynamic 
limit and in a transverse vector-field^^iS. This quan- 
tity vanishes for A'^ — > 00 in one dimension, but it can 
be nonzero for finite systems^. In Ref. 0, we have per- 
formed a study of the finite-size scaling of both quanti- 
ties for the XX Z chain, showing that « D^^ already 
holds at sufficiently high temperatures. At low tempera- 
tures and zero magnetic field, is always exponentially 
suppressed for even N due to finite-size gaps; thus a fi- 
nite value of Ds{T = 0) can only be found for 00. 
On the contrary, since D^^ w {■k/N){—T) at low temper- 
atures, correctly results in a finite value at T = 
in the massless regime. Depending on the context, one 
should carefully check which of these two quantities ex- 
hibits the more reliable finite-size behavior, and in fact, 
in the present case of finite magnetic fields we will argue 
in Sec. II VI that Dl should preferably be used. For a more 
detailed discussion of the relation between Dl and Bg^, 
we refer the reader to Ref. and references therein. 

In our numerical analysis, we will evaluate -Dth, £'s, 
and I?tii,s while the coefficients Dij from Eq. (|13|) can 
be derived if desired as they are linear combinations of 
Ah, A, and Dth.s- 

Dn = A, (18) 
Ai = Dth.s-hDs, (19) 
D22 = Ah-2/3/iAh,s + /3/i'A. (20) 

The XXZ model is integrable and solvable via the 
Bethe ansatz. Therefore one expects all quantities in 
Eqs. H13|) and (|14|) to be accessible by analytical tech- 
niques. Yet, for the spin Drude weight Ds{h = 0,r) at 
zero magnetic field, partly contradicting results can be 
found in the literature regarding both its temperature 
dependence and the question whether it is finite or not 
for the Heisenb erg chain (A = 1 ) in th e thermodynamic 
limit. See Refs. | 
cnces therein. 



and further refer- 



Ill. MEAN-FIELD APPROXIMATION 

We now discuss a Hartree-Fock type of approximation 
to the Hamiltonian Eq. which we use to compute 
the Drude weights Dij . The spin operators ' arc first 
mapped onto spinless fermions via the Jordan- Wigner 
transformation^^ 



(21) 
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Here, cl^-* destroys (creates) a fermion on site I. The string 

operator (f>; reads = X]i=i "-i with = c|cj. Next, 
the interaction term An;n;+i appearing in the fermionic 
representation is treated by a Hartree-Fock decomposi- 
tion leading to an effective mean-field Hamiltonian 



MF 



(22) 



with the mean-field dispersion 

efc - - J{(1 + 2Aa) cos(fc) + h~ 2A{n - 1/2)} . (23) 

The quantities to be determined self-consistently are a = 
(c[_(_]^C;) and n = {c\ci) where the latter is related to the 



average local magnetization m via m = {S[ 
The Drude weights can then be obtained from 

Dn = {■Kp/N){:il)- 
D21 = (7r/3/iV)(j2ji); 
D22 = {nl3yN){ji). 



- 1/2. 

(24) 
(25) 
(26) 



The current operators read 



Ji 



k 



J2 



with Vk = dck/dk. 

While this approach is exact for A = 0, fair results 
for Kti^{h = 0,T) are even obtained for < A < 1; 
see Refs. H andU From Eqs. ^ to the leading 

contribution at low temperatures can be derived. 

We start with the free fermion case A = 0, for which 
we find at the saturation field hc2 



with 



An 
A21 
A22 



4^ ] for h^h. 

All 



|(1-a/2)C(1/2), 



c2 



(27) 



3 

4 

15 /tt 
16 



|(2-^/2)C(3/2), 



(4-^/2)C(5/2); 



(28) 



C{x) being the Riemann-Zeta function. Note that the 
spin Drude weight at T = is finite for < ft. < hc2 
and vanishes for h > hc2- At low temperatures and 
for h = hc2, we find DnlT) = Ah^/T and a di- 
vergence of the pure thermal Drude weight Dth with 
-Dth ~ ftc2^ii^~^^^ to leading order in temperature, 
which follows from Eqs. (|19|) and H2U|I . We mention that 
the result D22 oc T'^/^ at the critical field was also found 
within a continuum theory suggested to describe trans- 
port properties of two-leg spin ladders>i£ 

In the intermediate regime, i.e., the gapless state (ii) 
[see Table n. 



A'th(ft,T) 



-v{h)T; v{h) ^J^/l-h? (29) 



holds at low temperatures, because the dispersion is lin- 
ear in the vicinity of the Fermi level for fcp 7^ 0, tt. Note 
that A'th(/i,T) w D22{h,T) for small T in this regime. 
Equation results in A'th = tt^ JT/3 for h = 0, which 
is, e.g., known from Ref. 16;. 

For \h\ > \hc2\ = 1, both D22 and the second term in 
Eq. Iini), i.e., Dj^/iTDn), are given by 



D22^DiJiTDii) = J-G 



-G/T 



T 



(30) 



to leading order in temperature and for T G, where 
G/J = \h\ — 1 is the gap. This implies that Kth{h,T) 
is strongly suppressed at low temperatures due to the 
cancellation of the contributions to Kti^{h, T) in Eq. H13|l . 
In fact, such cancellation occurs in the next-to- leading 
order in T as well. One can further show, taking into 
account the first non-vanishing contribution to /•sTth in 
Eq. lO, that 



3 

4 



(31) 



describes the low-temperature behavior of the thermal 
Drude weight above hc2- In Ref. 0, it has been ar- 
gued that 1522 oc cxp(-G/r)/v^ is a generic feature 
of gapped systems with a finite thermal Drude weight. 

We further point out that the ratio of the thermal 
Drude weight i^th and the spin Drude weight Dg ful- 
fills a Wiedemann-Franz type of relation in the low- 
temperature limit in all three cases, i.e., in the massless 
and the fully polarized state as well as for h = hc2- 



Ds 



LoT. 



(32) 



The constant Lq takes different values in the regimes (i) 
and (ii), but for the free- fermion case (and within mean- 
field theory as well) it is independent of the magnetic field 
in the massless and fully polarized state, respectively. 

Before turning to the mean-field theory for A > 0, 
let us briefly discuss which results can be expected from 
conformal field theory for the massless state. The ex- 
pressions for the spin and thermal Drude weight D22 
and Dn have the same structure as at zero magnetic 
field, i.e., D22{h,T) = {n^ /3)v{h, A)T and Dii{h,T) = 
K{h, A) v{h, A) with field-dependent velocity v and Lut- 
tinger parameter K (see, e.g., Refs. Hand 0)- This im- 
plies that the constant Lq appearing in Eq. H32|l is field 
dependent in the massless regime (see Ref.il6ifor h = 0): 



3K{h,A) 



(33) 



Furthermore, D21 vanishes in the continuum limit due to 
particle-hole symmetry. While a finite magnetic field ini- 
tially breaks this symmetry for the original bosonic fields, 
the original form of the Luttinger-liquid Hamiltonian is 
restored by introducing a shifted bosonic field^S,. This 
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FIG. 1: Thermal Drude weight Kth{h, T) of the XX Z chain 
for A = 0.1: comparison of mean-field theory (MF) and exact 
diagonalization (ED). The thermal Drude weight Kth{h, T) is 
shown for h, = 0,0.5, 1.1, 1.5. Increasing field is indicated by 
the arrow. Thick lines denote results from the Hartree-Fock 
approximation; thin lines: ED for A'^ = 18. Deviations at 
low temperatures for h = and h = 0.5 are due to finite-size 
effects in the ED results. 

has an interesting consequence for the low temperature 
behavior of the pure thermal Drude weight Dth- Namely, 
by solving Eqs. H19fl and H2U|I for Dth, one obtains 

Ah = D22 + -j^ = y t'T + , (34) 

which implies that Dth diverges at low temperatures with 
T^^ in the massless regime, consistent with results of 
Ref.lli 

Additionally, one obtains ifth in the massless regime 
and in the low-temperature limit 

ifth(/i,T)«D22 = y«(/i,A)T. (35) 

Both parameters, i.e., K = K{h,A) and v = v{h,A), 
can be computed exactly by solving the Bethc-ansatz 
equations^. The velocity v — v{h) has been calculated 
for A = 1 in Ref. Further numerical values for these 
parameters can be found in, e.g., Ref. lill 

Let us next dicuss the results from the mean-field ap- 
proximation (MF) for A > 0. Figure U] shows Kth{h,T) 
for A = 0.1 and h = 0,0.5,1.1,1.5 (thick lines). The 
main features are: (i) a suppression of the thermal Drude 
weight by the magnetic field; (ii) a shift of the maximum 
to higher temperatures for h > 0.5 compared to h ~ 0; 
(iii) a change in the low-temperature behavior which will 
be discussed in more detail below in this section. 

For comparison, the results from exact diagonalization 
(ED) for = 18 sites are included in Fig. ^ (thin lines) 
and we find that the agreement is very good. Devia- 
tions at low temperatures for h = and h = 0.5 are 
due to finite-size effects, i.e., the ED results are not yet 
converged to the thermodynamic limit. For larger fields 



h > hc2 = 1.1, deviations between ED and MF are neg- 
ligibly small. 

From Eq. (|23|) . we can derive the critical field hc2 
within the Hartree-Fock approximation. At T = and 
h = hc2, the ground state is the fully polarized state with 
n = (c|cj) = 1, i.e., the parameter a from Eq. H23|) van- 
ishes. Consequently, we find hc2 = 1 + A in accordance 
with the exact result. Indeed, the low-energy theories 
along the line h = 1 + A and for A = are equiva- 
lent in the sense that they are characterized by the same 
Luttinger parameteri^SiiS. Within bosonization, the line 
h = hc2 is particular since the velocity of the elementary 
excitations vanishes here. 

Regarding the low-temperature behavior of the ther- 
mal Drude weight we can then conjecture that it is given 
by Eqs. lf77|l and (^5)) for h = hc2, independently of 
A. We will come back to this issue in Sec. IIVI where 
we discuss the results from exact diagonalization for 
A > 0. The case of A — —1 and h ~ 0, however, 
seems to be an exception as we have found indications 
for Kth{h = 0) oc T at low temperatures before^. Here, 
the existence of many low-lying excitations might com- 
plicate the situation. 

In the ferromagnetic state and for low temperatures, 
the parameter a from Eq. (|23|) is exponentially sup- 
pressed and the average local magnetization is m = 1/2. 
Thus, to leading order in T the low-temperature de- 
pendence of Kthih,T) is independent of A, similar to 
the case of /i = hc2, and the thermal Drude weight is 
exponentially suppressed Kth{h,T) oc T^/^e~'^^'^ with 
G ~ h ~ hc2 ■ 

In the gaplcss state our mean-field theory results con- 
firm that Kthih, T) = V{h, A) T for A > and low tem- 
peratures. However, the mean-field prefactor V{h^ A) 
will be renormalized if interactions are fully accounted 
for; see Eq. 

In summary, we have obtained the leading low- 
temperature contributions to Kth in the regimes (i) and 
(ii) of TableQlusing mean-field theory and conformal field 
theory. Mean-field theory provides a reasonable quanti- 
tative description of the transport coefficients for small 
A and h as well as for h > hc2- 



IV. EXACT DIAGONALIZATION 

In this section, we first present numerical results for the 
thermal Drude weight of the Heisenberg chain (A = 1). 
Second, the field dependence of Kth{h,T) for interme- 
diate temperatures T is analyzed. Next, Kthih, T) for 
h = hc2 is discussed for different choices of the anisotropy 
A > and finally, we make some remarks on the lower 
bound for the spin Drude weight Dn = Ds given in 
Eq. (|ll|l . While Ds{h,T) still eludes an exact analytical 
treatment for arbitrary temperatures, analytically exact 
results for Dth{h,T) and Dth,s{h,T) of the Heisenberg 
chain have very recently been reported in Ref. ITqI 

Let us first address a technical issue, namely the ap- 
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FIG. 2: Thermal Drude weight Kii^{h,T) of the Heisenberg 
chain (A = 1). Panel (a): Kthih,T) computed from Eq. l(TI|l 
using Ds — (solid lines) and Dg ~ Dl' (dashed lines), 
both for h = 0.5 and iV = 12, 14, 16, 18, 20. Arrows indicate 
increasing system size. The dotted line is the result for /i = 
and N ^ oo from Ref. H^. Inset: comparison of Dl{h,T) 
[solid lines] and Di'{h,T) [dashed lines]; h = 0.5, A = 1. 
Panel (b): Kth{h,T) for h = 1.5,2,2.5,3 and iV = 20 [thick 
lines; Ds{h,T) = Di{h,T)]. The curve for N = 18,h = 1.5 is 
included (thin solid line). 



propriatc choice for Ds{h,T) in Eq. (|14|l . For the case 
of zero magnetic field, we know from our previous study 
Ref. that Dl{h,T) and D^J{h,T) exhibit a different 
finite-size behavior at /i = 0. This is similar to the situ- 
ation at finite fields. The inset of Fig. [2Ja) shows both 
Dl{h,T) and D^/{h,T) for A = 1 and /i = 0.5, and we 
see that first, D^^ {h, T) is well converged at low temper- 
atures; and second, a large difference between {h, T) 
and D^^{h,T) is visible at low temperatures. The ther- 
mal Drude weight Kth{h,T), resulting from either in- 
serting Dl{h,T) or Dl\h,T) in Eq. is shown in 
Fig. |2Ia). We have decided to use Dl in the numeri- 
cal study for consistency reasons, since then, all Drude 
weights entering in Eq. (|14|l have a similar finite-size de- 
pendence at low temperatures, characterized by the ex- 
ponential suppression at low temperatures due to the 
finite-size gap. On the contrary, using D^^{h,T) leads 
to an artificial double peak structure in A'th(/i,T); seen 
in Fig. Ha). 

We have checked that a similar scenario arises for 
A = for finite systems. However, for this case the 
Drude weight can be computed exactly in the thermo- 
dynamic limit and we find that one of the two maxima 
disappears. Thus we expect an analogous behavior for 
A > 0, supporting the choice of instead of D^^ . 

Further numerical results for Kt\^(h^T) of the Heisen- 
berg chain are provided in Fig. [^b) for h > 1.5. The 
main features of the thermal Drude weight can be sum- 



T/J=0.5, N=18 




^02 

FIG. 3: Main panel: field dependence of the thermal Drude 
weight Ktb{h,T) for A = 0.5,1,2 and T/J = 0.5 (ED for 
N — 18 sites). Inset: field dependence of Dti^{h,T) for the 
same parameter sets as in the main panel. 

marized as follows: (i) for < /i < hc2, finite-size effects 
arc small for T/J > 0.4 [see Fig.EIa)]; (u) for h > hc2, 
finite-size effects are negligible; (iii) the position of the 
maximum depends on the magnetic field; (iv) Kthih,T) 
is strongly suppressed as the magnetic field is increased. 

As both Dt\i{h, T) and Z)th,s(/i, T) converge rapidly to 
the thermodynamic limit at high temperatures, the small 
finite-size effects observed for T/J > 0.4 [see Fig. EJa)] 
arc due to Ds{h,T) [see the inset of Fig. I^^a)]. At low 
T, Kth{h,T) increases with system size N while it de- 
creases with growing N at high temperatures. The van- 
ishing of pronounced finite-size effects upon approach- 
ing the line h = hc2 from below can be ascribed to the 
fact that a description in terms of free fermions with 
parameters independent of A is valid here, as was al- 
ready evidenced in the previous section. For the ferro- 
magnetic state {h > hc2), the curves shown in Fig. |2Ib) 
for = 20 are indistinguishable from the corresponding 
ones for A^ = 18 (not included in the figure) within the 
line width. 

Regarding the position of the maximum, there is evi- 
dence that it is first shifted to higher temperatures when 
the field is switched on as compared to the case of ft, = 0; 
see Fig. |2Ia). A precise determination of its position in 
the intermediate gapless phase is somewhat complicated 
as typically, the numerical data converge well down to 
roughly only the peak temperature. Still, there are indi- 
cations that at strong fields /i ~ 1, the maximum tends 
to be located at lower temperatures than for h = 0. 
This can be seen, for instance, in the case of ft = 1.5 
in Fig.|2Ib). In the polarized state, Kt\i{h,T) definitely 
peaks at larger temperatures than at vanishing field due 
to its exponential suppression at low temperatures. 

The decrease of Kt\i{h,T) as a fimction of increasing 
magnetic field as mentioned in the preceding discussion of 
the Heisenberg chain is also observed for other choices for 
the anisotropy A. This is demonstrated for A = 0.5, 1, 2 
at T/J = 0.5 in the main panel of Fig. O where Kt\i{h, T) 
is shown as a function of the magnetic field ft and plotted 
versus ft/ftc2. In contrast to Kthih,T), Dtii(h,T) grows 
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FIG. 4: Thermal Drude weight Kth{h,T) at the critical field 
fec2 = 1 + A for A = 0,0.1,0.5, 1,2. For A / 0, we show 
numerical results for A*' = 18 sites, while the curve for the free 
fermion case (thin solid line) is valid in the thermodynamic 
limit. 



with increasing magnetic field at intermediate tempera- 
tures. For illustration, Dth{h,T) is plotted versus h/hc2 
at T/J = 0.5 in the inset of Fig. |21for the same choice 
of parameters as in the main panel, ft exhibits a maxi- 
mum at large fields, which increases and its position ap- 
proaches h = hc2 when the temperature is lowered. Thus, 
indications of the transition to the ferromagnetic phase 
arc visible in DiY^{h,T), but not present in Ki\^{h,T). 
Note, however, that all three curves in the main panel of 
Fig. O almost pass through the same point for h w hc2- 
For T/J < 1, Dth{h,T) is enhanced by the magnetic 
field and we can therefore conclude that the decrease of 
Kth{h, T) as a function of magnetic field is due to a can- 
cellation of Dth{h, T) and the magnetothermal correction 
in Eq. (HH). 

Along the critical line h — hc2 = 1 + A, further ev- 
idence for universal low-temperature behavior can be 
found by ED. This can be seen in Fig.^ where we present 
A'th(/i,T) for A = 0.1,0.5, 1,2 and N = 18. The curve 
for A = is also included in the figure; this one, however, 
is exact in the thermodynamic limit. Below T/J « 0.25, 
the curves lie on top of each other. Small deviations at 
lowest temperatures visible in the plot can be ascribed 
to the presence of finite-size gaps. This supports our 
conclusion from Sec. [ffl] that Eqs. (|?7|l and if^ hold for 
arbitrary A > and further numerical data (not included 
in the figure) show that it is also correct for —1 < A < 0. 

Finally, let us turn to the inequality Eq. for the 
spin Drude weight Dii(h,T) = Ds{h,T) introduced in 
Sec. |n] Here, we want to discuss to which extent the 
inequality Eq. (|ll|l is exhausted by jth at finite magnetic 
fields and finite temperatures. An analogous analysis in 
the limit of /3 = can be found in Ref. |22- To this end 
we compare Dl {h,T) and 
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in Fig. El Note that first, the relation Ds{h,T) > 
DsuhihyT) is equivalent to the positivity of the ther- 



FIG. 5: Comparison of the spin Drude weig ht Dl{h,T) (thick 
lines) and its lower bound Dsuh (thin lines; see text for further 
details). The figure shows results for A = 1, = 20, and 



mal Drude weight Kt^{h,T) > 0. Second, D^{h,T) w 
Dsuh{h,T) implies a very small thermal Drude weight 
and thus, the comparison provided in Fig. |5l also reveals 
the relative size of the two contributions to Kt\i{h,T) 
in Eq. H14() . namely Dth{h,T) and the magnetothermal 
correction D^^ ^{h, T)/[TDl{h, T)]. In Fig.gl results are 
shown for A = 1, TV = 20 sites, and h = 0.5,2,2.5. For 
the sake of clarity, data for smaller system sizes are not 
included in the figure. Differences between the curves 
for = 18 and = 20 are anyway only pronounced 
for temperatures T/J< 0.1 and become smaller as the 
magnetic field h increases. 

Figure |S1 allows for three major observations: (i) 
Dl[h,T) w Dsuh{h,T) at low temperatures and for all 
cases shown in the figure; (ii) Dsuhih,T) approximates 
Dl[h,T) the better the larger the magnetic field is; (iii) 
significant deviations are present for high temperatures 
implying that for a quantitative description of Ds{h,T) 
using Eq. (|ll|l . more conserved quantities need to be con- 
sidered in Eq. (|ll|l . 

Our comparison provides, at least for finite system 
sizes, a quantitative measure of the temperature range 
where « Dgub- Point (i) indicates that analytical 
approaches can make use of Dsuhih,T) for a quantita- 
tive description of Ds{h, T) at low temperatures as it has 
been done by Fujimoto and Kawakami within a contin- 
uum theory in Ref. l28l The quantities that appear on 
the right hand side of Eq. (|36|l are less involved than 
Eqs. H15|) and (|16|l for Ds{h,T), as the former arc static 
correlators. Furthermore, for finite magnetic fields, we 
suggest to compute Ds{h,T) analytically from Eq. Hll|l . 
taking into account some of the conserved quantities Qm, 
which are in principle known (see, e.g., Ref. l43l) . Such 
a procedure is applicable to /i ^ and might circum- 
vent the ambiguities in the result s encounte red in recent 
computations of D^ih = 0) fRefs. I28l36i38ll4l . The lat- 
ter have used Eq. Ijltill directly or Kohn's formula™*^. 
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equivalently. Regarding the relative size of Dth{h, T) and 
the magnetothermal correction, we see that the latter be- 
comes more relevant the larger the magnetic field is which 
leads to the strong suppression of Kth{h, T). This is con- 
sistent with results of the previous sections of this paper. 

V. CONCLUSIONS 

In this paper we have studied the thermal Drudc 
weight of the XX Z model with exchange anisotropy 
A > in finite magnetic fields using mean-field theory 
and exact diagonalization. Magnetothermal effects have 
been taken into account and the condition of zero mag- 
netization current flow has been applied. Let us now 
summarize the main findings and relate them to experi- 
ments. 

We have discussed the low-temperature limit of the 
thermal Drudc weight Kth{h,T) and we have given ar- 
guments that it changes from an algebraic behavior for 
< A < < hc2 to an exponentially activated be- 
havior in the polarized state for h > hc2- In addi- 
tion, the leading term at low temperatures along the 
critical line h = hc2,^ > —1, is universally given by 
Kth{h,T) = AT^/^, where the prcfactor A, given m 
Eqs. H27|l and H28|l . is independent of A. In the gapless 
phase, the leading contribution to Kt\i{h,T) is linear in 
the temperature with a field- and anisotropy dependent 
prefactor. In consequence, the thermal Drudc weight 
A'th(ft., T) can be expected to be proportional to the spe- 
cific heat in the gapless state in the low-temperature 
limit; where the velocity of elementary excitations is con- 
stant. 

Further, the Drude weight is suppressed by the mag- 
netic field, which can be ascribed to the increase of the 
magnetothermal correction relative to the pure thermal 
Drude weight Dti^{h,T). As a third result, the position 
of the maximum of i^th(^, T) depends non-monotonically 
on the magnetic field. While in the present paper, we 
have focused on the thermal Drudc weight A'th under 
the condition of zero spin-current flow, our analysis of 
the full transport matrix Eq. (j^J can easily be extended 
to a variety of other transport situations, which would 



be characterized by different combinations of the Drude 
weights. 

Turning now to experiments, we emphasize that for a 
description of realistic materials external scattering has 
to be accounted for. In a simple picture, one may expect 
the Drude peak to be broadened by external scattering 
mechanisms. The behavior in magnetic fields that one 
may observe in experiments will likely depend both on 
external scattering rates as well as on the thermal Drude 
weight. Nevertheless, one may speculate that qualita- 
tive trends of the field dependence of the thermal Drude 
weight are reflected in thermal transport experiments. 

In recent experiments on quasi one-dimensional mag- 
netic materialsiiSiSii^, the thermal conductivity has of- 
ten been found to be insensitive to the application of 
an external magnetic field. This is, however, explained 
by the large absolute value of the exchange coupling in 
these materials, being typically of the order of magnitude 
of 1000 K. It would therefore be desirable to perform 
measurements on materials with a moderately small ex- 
change coupling to check our results. Still, it might be 
very difficult to reach the saturation field hc2 in realistic 
materials, but at least the qualitative features like the 
suppression of the thermal conductivity or a shift of the 
maximum could be verified. In particular, experiments 
in relatively large magnetic fields may provide an indi- 
rect probe of spin currents in one-dimensional quantum 
magnets, which have so far not yet been observed directly 
in experiments. Conclusions about the low-temperature 
limit would require reliable methods to separate the mag- 
netic contribution from the phonon part, which is a chal- 
lenging task in the interpretation of experiments. Never- 
theless, we believe that further experiments on the field 
dependence of the thermal conductivity could hint at the 
nature and mechanisms of magnetic transport properties. 
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